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Abstract 

Using a realistic model of activity dependent dynamical synapses and 
a standard integrate and fire neuron model we study, both analytically 
and numerically, the conditions in which a postsynaptic neuron efficiently 
detects temporal coincidences of spikes arriving at certain frequency from 
A*' different afferents. We extend a previous work that only considers 
synaptic depression as the most important mechanism in the transmission 
of information through synapses, to a more general situation including 
also synaptic facilitation. Our study shows that: 1) facilitation enhances 
the detection of correlated signals arriving from a subset of presynaptic 
excitatory neurons, with different degrees of correlation among this subset, 
and 2) the presence of facilitation allows for a better detection of firing 
rate changes. Finally, we also observed that facilitation determines the 
existence of an optimal input frequency which allows the best performance 
for a wide (maximum) range of the neuron firing threshold. This optimal 
frequency can be controlled by means of facilitation parameters. 

1 Introduction 

Recently it has been reported that postsynaptic potentials recorded in cor- 
tical neurons present dynamical properties depending on the presynaptic ac- 
tivity l |Tsodyks and Markram, 1997| [Abbott et al., 1997| ). This behaviour can 
be understood by means of the action of different mechanisms taken place at 
the level of the synapses, as for instance, short-term depression and/or facil- 
itation. The first mechanism considers that the amount of available neuro- 
transmitters in the synaptic button is limited. Thus, the neuron needs some 
time to recover these synaptic resources in order to transmit the next incom- 
ing spike. As a consequence, the dynamics of the synapse is an activity- 
dependent mechanism producing a non-trivial effect in the postsynaptic re- 
sponse. This picture differs from the classical synaptic description which consid- 
ers the synaptic strengths as static identities with the only possible time mod- 
ification due to slow learning processes HHopfield, 1982| ). Moreover, it is well 
known that short-term depression plays an important role in several emerging 
phenomena in the brain such as selective attention ( Buia and Tiesinga, 2005] 



[McAdams and Maunsell, 1999| ), cortical gain control ( ^bo^^^^l997| ), and 
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complex switching behaviour between activity patterns in neural network mod- 
els dPantic et al., 2002) [Cortes et al., 2006| ). However, the complete study of 
other mechanism present in synaptic transmission on pyramidal neurons, as 
synaptic facilitation, that compete with depression still is lacking. Synaptic 
facilitation takes into account that the influx of calcium ions through voltage- 
sensitive channels favours the neurotransmitter vesicle depletion. This yields to 
relevant behaviour in synchrony and selective attentio n ||Buia and Tiesinga, 2005| ) 
and in detection of bursts of action potentials (AP) ( [Matvee'T^Ln^^^ingT^OOOl 
IDestexhe and Marder, 20041 ). Therefore, one would expect that synaptic facili- 
tation have a positive effect in the efficient transmission of temporal correlations 
between spikes trains arriving from different synapses. 

In this work, we used the phenomenological model of dynamic synapses, in- 
cluding both depressing and facilitating mechanisms, introduced in ( |Tsodyks et al., 1998| ) 
to study the cooperative effect of both in spike coincidence detection (CD) tasks. 
That is, we computed the regions, in the space of the relevant parameters, 
in which a postsynaptic neuron can efficiently detect temporal coincidences of 
spikes arriving from N different afferents. The aim is to determine the range 
of the parameters defining the dynamic of the synapses and neuron for which 
the performance of the neural system under study is improved. Our study 
shows that facilitation enhances the detection of correlated spikes and firing 
rate changes in situations for which the mechanism of depression alone does not 
perform well. These main results are robust and persist even when one decreases 
the degree of correlation between the afferents. Moreover, synaptic facilitation 
determines the existence of an optimal frequency which allows the best per- 
formance for a wide range of the neuron firing threshold. The location of this 
optimal frequency can also be controlled by means of facilitation control param- 
eters. This property can be important in neural media constituted by neurons 
presenting heterogeneity in the firing threshold ( |Azouz and Gray, 2000| ) in order 
to efficiently process information codified, for instance, at this frequency. 



2 The model 

We consider a postsynaptic neuron which receives signals from N presynaptic 
neurons through excitatory synapses. As a first approximation to model exper- 
imental data, we assume that the stimulus received by a particular neuron, as a 
consequence of the overall neural activity, is modeled by a spike train following 
a Poisson distribution with mean frequency / ||Tsodyks et al., 1998|). Accor ding 
to the phenomenological model presented in ( |Tso^y^^Ln^Ma^^n7l997| ), we 
consider that the state of the synapse i is governed by the system of equations 
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where Xi,yi,Zi are the fraction of neurotransmitters in a recovered, active and 
inactive state, respectively. Here, Tin and Tree are the inactivation and recov- 
ery time constants, respectively. Depressing synapses are obtained for U{t) = 
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Use constant, which represents the maximum amount of neurotransmitters 
which can be released (activated) after the arrive of each presynaptic spike. 
The delta functions appearing in |^ take into account that an AP arrives 
to the synapse at fixed time t — tsp. Typical values of these parameters 
in cortical depressing synapses are Tin — 3 ms, Tree = 800 ms, and Use — 
0.5 dTsodyks and Markram, 1997| ). 

The synaptic facilitation mechanism can be introduced assuming that U(t) 
has its own dynamics related with the release of calcium from intracellular stores 
and the influx of calcium from the extracellular medium each time an AP arrives. 
Here, we consider the dynamics proposed in ( |Tsodyks and Markram, 1997| ) which 
assume 

Uit)=uit)il~UsE) + UsE (2) 

with 

^ ^ + [/^^[i _ u(t)W - tsp). (3) 

at T fae 

Here, u(t) is a dynamical variable which takes into account the influx of cal- 
cium ions into the neuron near the synapse through voltage-sensitive ion chan- 
nels ( [Bertram et al., 1996| ). These ions usually can bind to some acceptor which 
gates and facilitates the release of neurotransmitters. A typical value for the 
facilitation time constant is Tfae = 530 ms ( [Markram et al., 1998| |. Then, the 
variable U{t) represents the fraction of neurotransmitters that are being ac- 
tivated, either by the arriving of a presynaptic spike (Use) and by means of 
facihtating mechanisms {u{t){l — Use))- 

One can think that the postsynaptic current generated in a particular synapse 
is proportional to the fraction of neurotransmitters which are in the active 
state, that is, li = Ase Hi, where Ase is the maximum postsynaptic cur- 
rent that can be generated ^. Hereafter, we will choose Ase ~ A2.5pA which 
is within the physiological range and gives an optimal system performance for 
Vth — 13mV, which is very near to the mean value threshold measured in some 
cortical areas ( [Azouz and Gray, 2000| ). Then, the total postsynaptic current 
generated by signals arriving from the N excitatory synapses can be written 
as Itotai — X^iLi ^i- This current generates a postsynaptic membrane potential 
which we modeled using an integration-and-fire (IF) neuron model, that is 

= - V + Rinl total, (4) 

where i?i„ = 0.1 Gfi and = 15 ms are, respectively, the input resistance and 
the membrane time constant. These typical values has been taken also from 
pyramidal cells ( [Tsodyks and Markram, 1997| ). The IF neuron model assumes 
that, once the membrane potential reaches a certain threshold Vth above the 
resting potential Vrest = 0, an AP is generated and V{t) is reset to zero. In 
addition, we assume the existence of a refractory period of r^e/ = 5 ms during 
which V{t) remains to zero after the generation of each postsynaptic AP. 

^Note that the synaptic conductance rather than the synaptic current depends on AsE'Vit), 
however our assumption for the current is a good approximation when the membrane potential 
V{t) is below the firing threshold Vth ^^nd Tm S> Ti„, so that V(t) remains constant during 
the time in which the synaptic conductance varies. 
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3 Detection of strongly correlated signals 



First, we have studied the postsynaptic response of a neuron receiving input 
signals from N — 1000 excitatory synapses, with a subset of AI — 200 synapses 
stimulated by identical spike trains. These strongly correlated afferents fire 
spikes simultaneously, so we can consider them as a signal term. The remaining 
N — M synapses receive uncorrelated spike trains which constitute, therefore, 
a noisy background of activity superimposed to the signal. We have inves- 
tigated, both analytically and numerically, spike coincidence detection (CD) 
experiments. Our interest is to determine the values of the parameters charac- 
terizing the neuron and synapse models, for which the postsynaptic neuron can 
detect the embedding signal (that is, its response is strongly correlated with 
the input signal) . A typical CD experiment is showed in figure ^ The figure 
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Figure 1: Response of a postsynaptic neuron receiving Poisson spike train 
(Top panel) at frequency of 10 Hz from N = 1000 afferents through dynamic 
synapses. Left and right panels corresponds to the case of depressing and fa- 
cilitating synapses, respectively. In these simulation. Use takes the values 0.5 
(middle panel) and 0.01 (bottom panel), respectively, and the threshold is fixed 
at 13mV. The figure shows that facilitation enhances CD tasks for relatively 
low values of Use- 

clearly shows the effect of including facilitation compared with the situation in 
which only depression is considered. For high values of the parameter Use, 
the system presents a good performance in the CD of the incoming signals, in 
both cases. However, for small Use the detection of signal is improved in the 
presence of facilitating mechanisms. In fact, when Use takes low values, the 
contribution of depression to U (t) (which gives the strength of the synapse) be- 
comes irrelevant. Facilitation, however, still contributes to maintain U{t) highly 
enough to allow a good performance on the CD task. For even smaller values of 
Use, for instance lower than 0.007, the depressing-facilitating mechanism also 
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fails. In general, given a particular value of r/ac there exists a critical value 

of Use below which the model does not perform well, even inchiding facilitat- 
ing mechanisms. On the contrary, for higher values of Use and in the region 
in which the model has a good performance, the simulations do not show any 
differences when one includes facilitating mechanisms. The reason is that the 
facilitation term becomes irrelevant for high values oi Use and depression is 
the only mechanism contributing to the dynamics. However, still there is a re- 
gion in the space of relevant parameters -for intermediate values of Use~ where 
facilitating mechanisms allows for a better performance in the CD tasks. 
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Figure 2: Coincidence detection maps for a system with facilitating synapses. 
The values of Use were 0.002 (top), 0.05 (middle) and 0.5 (bottom). One can 
see that the effect of increasing Use is the spreading of the region of good 
CD (light zone) to the right. Simulations (left) confirm the analytical results 
(center). The right panels represent the same CD regions with only depressing 
synapses. In all cases r/ac was 530 ms. 

For a more general evaluation and quantification of the role of the facilitating 
mechanism, we computed the fraction of errors that occur in the detection of 
the presynaptic signal by the postsynaptic neuron as a function of the incoming 
frequency / and the neuron threshold Vth- These CD error maps give us a better 
perspective of the regions where one has a good performance in the space of the 
relevant parameters. Thus, for each pair (/, Vth), we computed in the stationary 
regime 1) the number of coincidence-input-events in the subset of M coincident 
afferents, namely Ninputsi 2) the number of output spikes in the postsynaptic 
neuron occurring immediately within a time-window of A = 5 ms after the 
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coincidence- input-events, Nhus, 3) the number of output spikes which are not 
hits, Nfaisesj and 4) the number of coincidence-input-events which did not result 
in output spikes N failures within the time window A ( |Pantic et al., 2003[ ). The 
fraction of errors is then defined as 

rvy 7-, ^ f ailures ~t~ -^falses /^\ 

VoError = — — . (5) 

inputs 

Analytical expressions for the quantities appearing in jSJ have been obtained 




Figure 3: The error function (as defined in the text) for a system with facil- 
itating synapses. From top to bottom and for fixed Use — 0.05, we consider 
'''fac = 50, 530, 1500 ms, respectively. The figure shows that the numerically 
computed error function from simulations (left panels) fits the analytically eval- 
uated (see appendix) error function (right panels). The figure also shows the 
positive effect of considering stronger facilitation for a better detection of the 
embedded signal. This is shown by the presence of a larger CD region with a 
small value for the error function (white areas). 

by integration of the model equations H4I1() and their derivation is explained in 
the appendix. 

We have computed both theoretical and numerical CD error functions. The 
results are showed in figure |21 The light area corresponds to regions where the 
postsynaptic neuron is able to efficiently detect the coincidence-input-events 
and to generate a postsynaptic response strongly correlated with the embedded 
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signal. On the other hand, dark areas are regions with a high percentage of 
errors. These errors can be produced, for instance, when N failures is large, which 
occurs for Vtu very large (grey areas), or when Nfaises increases, normally for 
small Vth in such a way that any current can produce a false event (black areas) . 
The figure also shows for fixed Tfacj the influence of Use in the behaviour of the 
system. Thus, when its value increases (from top to bottom) the width of the 
light area enlarges and spreads to the right, allowing a better CD for regions with 
high thresholds. The left panels correspond to numerical simulations whereas 
the central panels are the same error function evaluated using the analytical 
formulas, derived in the appendix, into equation ©. The figure shows the good 
agreement between theory and simulations. In the right panels we computed the 
same regions but considering only the mechanism of synaptic depression. One 
observes that for only depression and a limited amount of neurotransmitters 
{Use < 0.5), the region for good CD is narrower. In this situation, a large 
region of good detection is obtained only for Use near to one. 

Then, we concluded that for the same value of the amount of activated neu- 
rotransmitters the overall performance of the system is better with facilitation 
than if one only consider depressing mechanisms. This conclusion can be also 
observed when one fixes Use and varies the facilitation time constant Tfac, as it 
is shown in figure A large value for Tfac means an increase in the duration of 
the facilitating effect. As a consequence, the region for good detection enlarges 
compared with the situation of only depression, in special when the fraction of 
available resources is not too high. 

A detailed observation of figures 12 and |21 also shows the existence of a cer- 
tain frequency which allows a good performance for a wide (maximum) range 
of values of Vth (see for instance, middle panels in figure |21 that shows a good 
performance in detecting signal frequencies around 10 Hz for a threshold rang- 
ing from 8 to 18 mV). This optimal frequency decreases as Use goes to higher 
values, and becomes zero when only depression mechanism is relevant. There- 
fore, the presence of facilitation reveals the appearance of an optimal frequency 
and allows to control it by tuning the facilitation parameters. This result could 
be important to understand how real neural systems -where different types of 
neurons may have non-identical firing thresholds- can self-organize to efficiently 
detect and process correlated signals. 

Finally, to quantify the variation of the regions of good CD in the presence 
of facilitation and/or depression, we computed the maximum range of frequen- 
cies, A/, for which the postsynaptic neuron can detect signals with small error 
(less than 0.5). We choose a fixed threshold around Vth — HmV, and study the 
system behaviour for fixed Use and varying Tfac and vice versa. The results are 
presented in figure 01 The figure shows (left panel) that the range of good CD 
decreases with Use for the case of only depressing synapses, and even vanishes 
for Use < 0.05. However, if facilitating mechanisms are also present the sys- 
tem is able to recover the good performance by increasing the facilitation time 
constant. The figure also reveals (right panel) that facilitation always enlarges 
the maximum range of frequencies for any fixed value of Use (note that the 
depressing synapses limit is obtained for Tfac = 0). Thus, we conclude that for 
any values of these parameters the inclusion of facilitation in the dynamic of 
synapses improves the detection towards wider ranges of frequencies. Similar 
results are found for the maximum range of thresholds which allows good de- 
tection of a given frequency /, see El As a final conclusion, these results show 
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that facilitation expands the regions where the system can detect signals with 

small errors, both for a wide range of frequencies and neuron thresholds, and for 
any possible values of the relevant parameters defining the dynamic of synapses, 
namely Use and Tfac- 
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Figure 4: Variation of A/, as defined in the text, for different values of Use and 

Tfac- Left: A/ as a fimction of Tjac for three fixed values of Use- Right: A/ as 
a function of Use for three fixed values of Tfac- In all cases Vth was 13 mV. 
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Figure 5: Variation of AVth, as defined in the text, for fixed Use and variable 
Tfac (left) and vice versa (right). We can see that the case of Tjac = (depression 
limit) has always a worse behaviour in CD than the facilitation case. The 
frequency is set to 10 Hz. 



4 Effect of jitter 

Detection of coincident signals arriving from different presynaptic neurons have 
been treated in the previous section in an approximate way, that is, the embed- 
ding signal was constituted by fully correlated temporal events. However, in a 
real situation the incoming signals arriving to a neuron from different synapses 
are not totally correlated in time. Then and following our previous analysis. 



8 



if the signal term is formed by presynaptic spikes arriving from M afferents, 
the presynaptic neurons would never fire exactly at the same time to. On the 
contrary, it is more realistic to consider that the M presynaptic neurons will 
fire at random times U distributed around to following, for instance, a Gaussian 
distribution p{ti) with a certain deviation or jitter a. In this section, we consider 
the implications of this assumption to test the validity of the results previously 
obtained, and to investigate the effect of the jitter in the detection of signals 
that are not fully correlated. 




Figure 6: Excitatory postsynaptic current for a system receiving an single 
AP from M = 200 afferents. In each presynaptic neuron the AP occurs at 
different time ti which for all neurons is Gaussian distributed around to = 
300ms. The figure shows that the effect of jitter is the spreading of the current 
curve whereas an increment in the inactivation time constant causes longer right 
tails. Numerical results (lines) are in concordance with the analytical derivation 
of the current (dots) (see main text for an explanation). 

We start by computing the excitatory postsynaptic current generated in a 
synapse i due to a single presynaptic AP occurring at time ti, that is 

hit) = Ipeak exp[-(t - ti)/nn] t > ti (6) 

where Ipeak is the steady-state maximum current through a synapse obtained 
after stimulation with a periodic spike train (see the appendix for details) . Since 
ti is a Gaussian distributed stochastic variable with (ti) = to and standard 
deviation a, qi(t) = exp[— (t — ti)/Tin] (with t fixed) is also a random variable 
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with range [0, 1] and probability distribution given by 



2t- 



exp 
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where erfc(x) = 1 — erf(x) and erf(x) is the error function. One can easily 
compute the two firsts moments for P[qi(t)\ obtaining: 
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In the case of many afferents, the total current in the postsynaptic neuron is 
I(t) = Ipeak X^il*! ^^(Oj where 1 < < M is the fraction of the M afferents in 
which the AP has already generated a postsynaptic response at time t, and it is 
given by i'{t) = M J_^p{ti)dti. This number depends on time due to existence 
of the jitter that desynchronizes the arriving of the AP in all afferents. Then 
for t -^to i'{t) is small, but for t near to and large than tg, ^{t) is high and we 
can use the central limit theorem to obtain: 

I{t) ^ Ipeak^{t,to), (10) 

where ^{t, tg) is a Gaussian variable with mean and variance given by 

mto))^Ht){q{t)) (11) 



^ .Jl.{t)[{q^t))~{q{t)f] (12) 

with 1/(1) — ^ [erf {^-^) + l] • Hereafter we will use this analytical approach 

to compute CD maps with a jittered signal. 

Since u(t) needs to be high in order to use the central limit theorem, one 
expects that the theoretical current defined by equations H10I12|I will fit better 
the numerical results for t > to. This is shown in figure where the analytically 
computed current after the arriving of M jittered APs (green dots) is compared 
with the simulated current (red curve), for different values of the jitter a and 
different values of the inactivation time constant Ti„. The figure shows the 
good agreement between the theoretically and numerically computed currents. 
Moreover, one observes that, the effect of increasing the jitter is the temporal 
spreading of the current so that the signal influence occurs during a large period 
of time but with a smaller amplitude. This will cause a small decreasing in the 
capacity of the system to detect spikes. On the other hand, if we fix the jitter the 
effect of increasing Tin is the appearance of longer tails for t > to, which would 
be a desirable effect since the response to the next incoming AP will be higher. 
However, no changes are detected in the amplitude of the current when Tin is 
modified. Note that the effect of jitter does not depend on other parameters 
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driving the dynamics of synapses, as Use, Tree or Tfac which only affect to the 
amphtude Ipeak (see the appendix). Therefore, one should not expect a strong 
effect of jitter on the emergent properties due to facilitation and/or depression. 

In order to compute CD maps we have to compute the voltage generated by 
the jittered signal, so that we have to integrate the Langevin equation 

= -V + Rinlpeak^iit^to) (13) 
to 

Here the sum extends to a train of events, each one consisting of M jittered 
AP centered around a particular instant of time Iq in the event's train. In order 
to give a first approximation to the solution of this equation, the fluctuations 
are neglected. Therefore the factor ^{t, to) is now a Gaussian function of time, 
centered at t^ for each event in the train. Using standard methods and assum- 
ing a periodic train of events occurring at t^ — 0, 1//, 2//, . . . , one can easily 
integrate the equation (|13(l obtaining 



V{t) = exp(-t/T„0 
where 



-p^-^/^--") i-crpi-vko ^^^ 



(14) 



W{t) = ^ I eMt'/r„,)I{t')dt' (15) 



m JQ 



I{t)= Ipeak {^{t,0)) (16) 

which determines the evolution of the membrane potential. Simulation shows 
that this expression is also valid for Poisson event trains (data not shown) . Then, 
one can use (I14|l to evaluate the CD maps similarly to the case of ct = (non- 
jittered events). Indeed, as it is shown in the appendix, to do that is necessary 
the evaluation of the maximum value of V(t) generated by the signal term (V^i) 
during the signal event duration. In the practice, this can be analytically done 
only in the case of cr = 0. For a ^ 0, Vm must be numerically computed from 

The maps for the detection of jittered events are presented in figure An 
important conclusion is that the CD maps here are qualitatively the same as 
the maps obtained previously in the zero-jitter case. Increasing the value of 
the jitter yields to a decreasing of the area of good performance, as one could 
expect. However, this effect in the light zone is not too dramatic. Indeed, one 
observes that the shape of the regions does not change too much for two very 
different values of the jitter, namely a = 0.5 (top panels) and cr = 3 (bottom 
panels) (see figured). The figure also shows the good agreement between theory 
(right panels) and simulations (left panels). In addition, the jitter also causes a 
small delay on the reaching of the membrane threshold (cf. figure where one 
has the event at = 300 ms and the maximum of the generated current occurs 
at t — to + 5t). This fact turns into an increment in the number of failures and 
false hits. Then, the numerical counting of hits, failures and falses is affected 
because we consider a fixed detection temporal window of A = 5 ms. For cr = 3, 
we solved this problem increasing A to 7 ms and the result is showed in figure[3 
(bottom panels). Such effect must, however, be taken into account in situations 
with a high value of the jitter. 
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Finally, the analysis presented in this section concludes that the main results 
obtained in the previous section are robust for a more realistic treatment of the 
input presynaptic signals. 
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Figure 7: Coincidence detection maps for the case of arriving presynaptic 
signals with a certain value of the jitter. A comparison between situations 
with low ((7 = 0.5) and high (cr = 3) jitter is showed (top and bottom panels 
respectively). One can see that the effect of increasing the jitter is a small and 
non-relevant decreasing of the good CD region. These results have been found 
both numerically (left) and analytically (right). The values for the facilitation 
parameters were Tfac = 1000 ms and Use = 0.1. 



5 Detection of presynaptic firing rate changes 

In the previous study we have considered the overall firing rate as a fixed 
parameter. This assumption is not realistic and, more interesting, is to con- 
sider the presynaptic firing rate as a dynamic variable as it happens in real 
neuronal tissue. The rate changes during normal functioning of neural sys- 
tems in the presynaptic current leads to a transient behaviour in the exci- 
tatory postsynaptic potential (EPSP) which could cause a burst or an AP 
in the postsynaptic neuron dTsodyks and Markram, 1997| [Pantic et al., 2003| 
[Abbott et al., 1997J . The question that arises is if the postsynaptic neuron 
is able to detect synchronous changes (increases) in the afferent firing rates. 
This property have been found only for depressing synapses and not for static 
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synapses ( |Pantic et al., 2003) . Another question is if synaptic facilitation could 
have some positive effect in the detection of these rates changes by the postsy- 
naptic neuron. In this section we try to answer this question by studying the 
effect of increasing facilitation in spite of depression in the optimal detection of 
rate changes in the presynaptic current. 

To start, we assume a population of TV = 1000 afferents firing uncorrelated 
Poisson spike trains with a certain frequency / into a postsynaptic neuron. This 
population changes its mean firing rate every 1000 ms. The figure |H1 shows a 
comparison in the output of the postsynaptic neuron for facilitating and depress- 
ing synapses. The threshold for firing was fixed in Vth = 17 ms and Use = 0.1. 
Simulations show that facilitating synapses (r/ac = 500 ms) allow for a better 
detection of rate changes, and over a large range of frequencies, than depressing 
synapses. In general, the regions in which depressing and facilitating synapses 
perform well can vary. Thus, there are particular situations where facilitation 
is needed to detect presynaptic rate changes. 

A simple theoretical approach can help us to find the regions of good firing 
rate changes detections. To obtain such transient behaviour which allows to 
a rate-change detection, the threshold of the postsynaptic neuron must satisfy 
C/2i^(/i) > Vth > Cf2iLj{f2), where /i is the initial rate, /2 is the firing rate 
after the change, C = RinNrin and is the stationary postsynaptic current 
strength for a given frequency. From the system of equations and following 
a reasoning similar to the strategy used in the appendix, one can easily obtain 

= Tj- (17) 

where Uoo is the steady state valued of U{t) (see the appendix). If now we fix 
the frequency step Sf = f2 — fi, the resulting expressions will only depend on 
/i. Since for large enough frequencies Cf2LLi{fi) is a decreasing function of /i 
and C f2Lo{f2) is an increasing function of /2 (and therefore of /i), these two 
tendencies will converge for some /i. This leads to a close area of good rate 
change detection between the two curves. 



6 Discussion 

In the last years, there was an increasing interest in the study of the computa- 
tional functionality of synaptic activity dependent processes, as synaptic facilita- 
tion and depression, in real systems and neural networks models ( Abbott a nd Regehr, 20041 



[Destexhe and Marder, 2004, ,Pantic et al., 2002, .Pantic et al., 2003J . In partic- 
ular, recently it has been reported how these processes affect the stability of 
the attractors of the dynamics and, as consequence, it emerges an oscillatory 
phase where the activity of the system is continously jumping among the attrac- 
tors dPantic et al., 2002| [Cortes et al., 20061 [Torres et al., 2006| |. This jumping 
behavior could explain the emergence of voltage transitions between up and 
down states observed in cortical areas ( [Holcman and Tsodyks, 2006) . 

In this paper, we have presented a detailed study of how the competition 
between synaptic facilitation and depression affects the neural detection of an 
embedded signal in a background of uncorrelated noise. Our study shows that 
the inclusion of the facilitation mechanisms enhances the performance of cor- 
tical neural systems. In particular, the transmission of information codified 
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Figure 8: Detection of firing rate changes with depressing and facihtating 
synapses. The top panel shows the mean firing rate of the N — 1000 presy- 
naptic neurons as a function of time. Middle and bottom panels shows the 
response of the postsynaptic membrane potential for facilitating and depress- 
ing cases, respectively. In these simulations parameters were Use = 0.1 and 
Tfac = 500ms(0 ms) for the facilitating (depressing) case, respectively. Detec- 
tion of variations onto lower frequencies are not possible with these synaptic 
mechanisms. 

in spike trains through the synapses is better and the detection of firing rate 
changes is also improved compared with the case of only depression. Thus, 
contrary to which it happens with only depression, the presence of facilitation 
makes not necessary to have a high value for the maximum amount of active 
neurotransmitters to efficiently detect correlated signals. This would lead us 
to think that facilitation has a crucial role in the processing of information 
through synapses even when the neuron does not have enough synaptic active 
neurotransmitters. We have also seen that, although important, it is not crucial 
to have a strong correlation between different presynaptic afferent to have a 
good detection of signals. Facilitation also determines the existence of an opti- 
mal frequency which allows good performance for a wide range of neuron firing 
thresholds. This result could be important to understand how real neural sys- 
tems -where different types of neurons with non-identical firing thresholds are 
connected in a complex way- can self-organize to efficiently detect and process 
relevant information ( |Azouz and Gray, 2000| ) . Our analysis has also shown that 
the main conclusions are also valid for realistic jittered signals. 
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Appendix: Analytical derivation of the error func- 
tion 

In this section we derived analytical expressions for the functions appearing in 
the definition of the error function ((SJ used to obtain theoretically the regions 
for good spike coincidence detection in the (/, Vth) parameter space. 

First, we assume that the total presynaptic current can be splitted in two 
terms: a signal term containing the correlated embedded signal and a noise term 
formed by the background of uncorrelated spikes. 



Noise contribution 

To take into account the noise generated by — M uncorrelated spikes trains, 
we assume that the current at time t = t* +t generated by a single spike arriving 
to the synapse i at time t* is given by 

/i(r, t*) = Ipeak exp(-T/Ti„) (18) 

where Ipeak represents the averaged stationary EPSC amplitude obtained after 
stimulation with a periodic spike train, assumption that we also suppose valid 
for Poisson distributed spike train. After this consideration, one easily obtains 
from equations H1I8|I that 

r . [/oo(l-exp(-l//T^ec)) ..„^ 

Vafe = AsE- — 7z — TT-^ — r~rre — r 

1 - (1 - Uoo) exp{-l/ fTrec) 

with Uoo = Moo(l — Use) + UsE, where Uoo is the value of u{t) in the stationary 
state (t oo). For a periodic spike train, Uqo is given by 

„ -TT exp(-l//T/ae) . . 

1 - (1 - Use) exp(-l//r/ac) 

We can compute the mean noise contribution of the current and fluctuations 
using the standard expressions 



-^noise — (-^) ; 

From these definitions and using the central limit theorem we obtain 



(21) 



/„o,,e ^{N- M)AsEfr.nU^ ^ ,1 , (22) 

1 - (1 - Uoo) exp(-l//Trec) 

where we assumed that Tin ^ Tree- 

If we neglect fluctuations {(^i^^-^^ — 0), we can write Vnoise = Rinlnoise- 
Using this expression one can compute Nfaises taken into account that false 
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firing occurs when Vnoise > Vth so by a direct integration of equation (0J in a 
period of time T gives Nfaises ~ T/ {r^e/ - r,„ ln(l - Vt h/Vnoise)} l|Koch, 1999| ). 
Now using that / = Ninputs/T, we finaUy obtain as in ( |Pantic et al., 2003^ 

\n{l-Vth/Vno^se)) ' ' 

where 9{x) is the Heaviside step function, which takes into account that for 

Vnoise < Vth Nfaises = 0. 

To take into account fluctuations of Inoise one can use the so caUed hazard 
function approximation ( |Plesser and Gerstner, 2000| ) but it has been reported 
that it gives the same results than those obtained using the formula H23() for 
high frequencies and, on the contrary to the expression H23(l . it does not work 
properly for small frequencies ( |Pantic et al., 2003| ). Therefore, hereafter we will 
neglect fluctuations in Inoise and use l|23|l as an approximatively valid expression 
to analytically compute Nfaises- 



Signal contribution 

To analyse the signal contribution (arising from M coincident spikes) we used the 
same method developed in ( |Pantic et al., 2003| ) for the case of only depressing 
synapses. That is, assuming that V{0;t*) is the membrane potential at t — t* 
when M coincident spikes arrive, by direct integration of the equation Q the 
membrane potential at time t = t* + t is 

V{T;t*) = e^/^"- lv{0;t*) + ^^lIlMIl^[e"r _ (34) 

L Tma J 

where a = "^j." and Ipeak is given by (|19|l including all the effects due to 
synaptic depression and facilitation. If the next signal event (M coincident 
spikes) occurs aX t = t' one can obtain the following recurrence relation: 

V{Q-t') = e^*/^'" \v{Q-t*) + ^^luMIpS^^e'^^t _ (35) 



with At = t' — t* . which allows for computing the stationary value for the 
membrane potential at the exact time of the signal event arrival (see also 
| |Kistler and van Hemmen, 1999| )), that is: 

-At/T„. RinMIpeak 6°^* - 1 

We define Vm gnai ss the maximum of the membrane potential reached between 
the arrival of two consecutive signal events separated by a time At. This can 
be easily computed from equation (|24|l with 1/(0, t*) replaced Vst- 



signal 



T„t(l - exp(-l//r„)) 
Tin{l - exp(-l//ri„)) 

where we consider t = At )^ 1//. 



RinMIpeak (27) 
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The expression of Vsignai allows for evaluate the number of failures assuming 



that N 



failures 



inputs 



hits • 



Then, one obtains by direct integration of 



equation an using the same reasoning that for N false case that 



N 



failures 



inputs 



1 - 



0{Vn, 



signal 



Vth) 



ref 



, ln(l - {Vth - Vsignai) /Vnoise)] 



(28) 



where we have considered a hit event every time Vn, 

< Vth we will have N 



■Vsignai reach Vth- Note 



that from if Vnoise + Vsignai < Vth we win have N failures — ^inputs • 

Expression for N falsest N failures allows for theoretically compute the number 
of errors in the CD maps. 
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